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A new two-parameter family of quasi-exactly solvable quartic polynomial potentials V(x) — — x 4 + 
2iax 3 + (a 2 — 2b)x' 2 + 2i(ab — J)x is introduced. Until now, it was believed that the lowest-degree 
one-dimensional quasi-exactly solvable polynomial potential is sextic. This belief is based on the 
assumption that the Hamiltonian must be Hermitian. However, it has recently been discovered 
that there are huge classes of non-Hermitian, "PT-symmetric Hamiltonians whose spectra are real, 
discrete, and bounded below. Replacing Hermiticity by the weaker condition of VT symmetry allows 
for new kinds of quasi-exactly solvable theories. The spectra of this family of quartic potentials 
discussed here are also real, discrete, and bounded below, and the quasi-exact portion of the spectra 
consists of the lowest J eigenvalues. These eigenvalues are the roots of a Jth-degree polynomial. 



PACS number(s): 03.65.Sq, 02.70.Hm, 02.90. +p 

Quantum-mechanical potentials are said to be quasi- 
exactly solvable (QES) if a finite portion of the energy 
spectrum and associated eigenf unctions can be found ex- 
actly and in closed form QES potentials depend on 
a parameter J; for positive integer values of J one can 
find exactly the first J eigenvalues and eigenfunctions, 
typically of a given parity. QES systems can be classified 
using an algebraic approach in which the Hamiltonian is 
expressed in terms of the generators of a Lie algebra ^ |) . 
This approach generalizes the dynamical-symmetry anal- 
ysis of exactly solvable quantum-mechanical systems, 
whose entire spectrum may be found in closed form by 
algebraic means ||. 

An especially simple and well known example of a QES 
potential M is 



V(x) = x e - (47-1) 



(1) 



The Schrodinger equation, — ip"(x) + \V(x) — E]ijj(x) = 0, 
has J even-parity solutions of the form 



tp(x) 



-x 4 /4 



J -I 

E 

fe=0 



c k x 



21, 



(2) 



The coefficients c k for < k < J — 1 satisfy the recursion 
relation 

4( J - fc)c fe _! + Ec k + 2(k + l)(2fc + l)cfc+i = 0, (3) 

where we define c_i = cj = 0. The simultane- 
ous linear equations (||) have a nontrivial solution for 
Co, ci, cj-i if the determinant of the coefficients van- 
ishes. For each integer J this determinant is a polyno- 
mial of degree J in the variable E. The roots of this 
polynomial are all real and are the J quasi-exact energy 
eigenvalues of the potential (pi). 



The lowest-degree one-dimensional QES polynomial 
potential that is discussed in the literature is sextic. How- 
ever, in this paper we introduce an entirely new two- 
parameter class of QES quartic polynomial potentials. 
The spectra of this family of potentials are real, discrete, 
and bounded below. Like the eigenvalues of the potential 
(Q), the lowest J eigenvalues of these potentials are the 
roots of a polynomial of degree J. 

The potentials introduced here have not been dis- 
covered so far because they are associated with non- 
Hermitian Hamiltonians. Recently, it has been found 
that there are large classes of non-Hermitian Hamilto- 
nians whose spectra are real and bounded below H|J. 
Although they are non-Hermitian, these Hamiltonians 
exhibit the weaker symmetry of VT invariance. A class 
of these Hamiltonians, 



H=p z -{ix) N (N>2), 



(4) 



was studied in Ref. ||]. The special case N = 4 corre- 
sponds to the Hamiltonian 



H 



V 



(5) 



It is not at all obvious that this Hamiltonian has a pos- 
itive, real, discrete spectrum. To verify this property, 
we must continue analytically the Schrodinger equation 
eigenvalue problem associated with H in (^) from the 
conventional harmonic oscillator (N = 2) to the case 
N — 4. In doing so, the boundary conditions at |x| = oo 
rotate into the complex x plane. At N = 4 the boundary 
conditions on the wave function tp(x) read 



lim ip(x) = 0, 

\x\ — >oo 



(6) 



where the limit x — > oo is taken inside two wedges 
bounded by the Stokes' lines of the differential equation. 
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The right wedge is bounded by the Stokes' lines at 0° and 
—60° and the left wedge is bounded by the Stokes' lines 
at —120° and —180°. The leading asymptotic behavior 
of the wave function is given by 

ip(x) ~ e~ tx3/3 (|a|->oo). (7) 

It is easy to see that the asymptotic conditions in @ 
are satisfied by ip(x). A complete discussion of the an- 
alytic continuation of eigenvalue problems into the com- 
plex plane is given in Ref. |l0| . Note that for all values of 
N between 2 and 4, the Hamiltonian (|4|) is not symmetric 
under parity. This parity noninvariance persists even at 
N = 4; eigenfunctions ip(x) of (|^) are not symmetric (or 
antisymmetric) under the replacement x — > —x. 

In this paper we generalize the Hamiltonian (]|) to the 
two-parameter class 

H = p 2 - x i + 2ma; 3 + (a 2 - 2b)x 2 + 2i(ab - J)x, (8) 

where a and 6 are real and J is a positive integer. The 
wave function ip(x) satisfies the boundary conditions (|^) 
and the differential equation 

E^(x) = -ii"{x) + [ - x 4 + 2iax 3 + (a 2 - 2b)x 2 

+ 2i(ab- J)x]ip(x). (9) 

We obtain the QES portion of the spectrum of H in 
(0) as follows. We make the ansatz 

i>{x) = e -^ 3 / 3 -° x2 / 2 - <6:E Pj_i(x), (10) 

where 

J -2 

Pj-i(x) =x J - 1 + Y / CkX k (11) 

is a polynomial in x of degree J — 1. Substituting tp(x) 
in ([h]) into the differential equation (^) , dividing off the 
exponential in ([To|), and collecting powers of x, we obtain 
a polynomial in x of degree J — 1. Setting the coefficients 
of x k (1 < k < J — 1) to gives a system of J — 1 
simultaneous linear equations for the coefficients Ck (0 < 
k < J — 2). We solve these equations and substitute 
the values of Ck into the coefficient of x°. This gives a 
polynomial Qj(E) of degree J in the energy eigenvalue 
E. The coefficients of this polynomial are functions of 
the parameters a and b of the Hamiltonian H in (j^). 
The first five of these polynomials are 

Qi = E - b 2 - a, 

Q 2 = E 2 - (2b 2 + Aa)E + 6 4 + 4a6 2 - 46 + 3a 2 , 
Q 3 =E 3 - (3b 2 + 9a)E 2 + (3b 4 + I8ab 2 - 166 + 23a 2 )£ 

-b 6 - 9ab 4 + 16b 3 - 23a 2 b 2 + 48ab - 15a 3 - 16, 
Q 4 = E 4 - (4b 2 + 16a)E 3 + (6b 4 + 48a6 2 - 406 

+86a 2 )^ 2 + (-46 6 - 48a6 4 + 806 3 - 172a 2 6 2 



+320a6 - 176a 3 - 96)£ + 6 8 + 16a6 6 - 406 5 
+86a 2 6 4 - 320a6 3 + 176a 3 6 2 + 2406 2 
-568a 2 6 + 105a 4 + 384a, 
Q 5 =E 5 - (5b 2 + 2ba)E 4 + (106 4 + 100a6 2 - 806 
+230a 2 )£; 3 + (-106 6 - 150a6 4 + 2406 3 
-690a 2 6 2 + 1200a6 - 950a 3 - 336)£ 2 + (56 8 
+100a6 6 - 2406 5 + 690a 2 6 4 - 2400a6 3 
+ 1900a 3 6 2 + 16966 2 - 5488a 2 6 + 1689a 4 
+3360a)E - 6 10 - 25a6 8 + 806 7 - 230a 2 6 6 
+1200a6 5 - 950a 3 6 4 - 136 06 4 + 54 88a 2 6 3 
-1689a 4 6 2 - 8480a6 2 + 7440a 3 6 + 30726 
-945a 5 - 7632a 2 . (12) 

The roots of Qj(E) are the QES portion of the spectrum 
of H. 

The polynomials Qj(E) simplify dramatically if we 
substitute 

E = F + b 2 + Ja (13) 

and 

K = Ab + a 2 . (14) 

The new polynomials have the form 

Qi = F, 

Q 2 =F 2 - K, 

Q 3 = F 3 - AKF - 16, 

Qi = F 4 - 10KF 2 - 96F + 9K 2 , 

Q 5 = F 5 - 20KF 3 - 336F 2 + 64K 2 F + 768K, 

Qe = F 6 - 35KF 4 - 896F 3 + 2h9K 2 F 2 + 7040KF 

-22bK 3 + 25600, 
Q 7 = F 7 - 56KF 5 - 2016F 4 + 784K 2 F 3 + 3h7\2KF 2 

-23MK 3 F + 288000P - 55296-ftT 2 , 
Qs = F 8 - MKF 6 - 4032F 5 + 1974K 2 F 4 

+ 132480KF 3 - 12916K 3 F 2 + 1760256F 2 

-681408if 2 F + 11025/Y 4 - 6322176X. (15) 

The roots of these polynomials are all real so long as 
K > -Kcriticah where identical is a function of J. The 
first few values of ^critical are listed in Table |. At K = 
-^critical the lowest two eigenvalues become degenerate 
and when K < A" cr iti ca i some of the eigenvalues of the 
QES spectrum are complex. Thus, the QES spectrum is 
entirely real above a parabolic shaped region in the (a, 6) 
plane bounded by the curve a 2 + 46 = identical- 
Extensive numerical calculations lead us to believe that 
the non-QES spectrum is entirely real throughout the 
(a, b) plane and that when K > A cr j t i C ai the eigenvalues 
of the QES spectrum lie below the eigenvalues of the non- 
QES spectrum. However, as we enter the region K < 



2 



^critical some of the eigenvalues of the QES spectrum 
pair off and become complex. Other eigenvalues of the 
QES spectrum may cross above the eigenvalues of the 
non-QES spectrum. In Fig. [l]we illustrate the case J = 3 
and o = 0. Note that for b > | the QES eigenvalues are 
three lowest eigenvalues of the spectrum. When b goes 
below §, two of the QES eigenvalues become complex and 
the third moves into the midst of the non-QES spectrum. 

The standard way to understand QES theories is to 
demonstrate that the Hamiltonian can be expressed in 



LU 




- 



terms of generators of a Lie algebra. Following Turbiner 
Q , we use the generators of a finite dimensional represen- 
tation of the SL(2, Q) with spin J. The three generators 
have the form 



If we apply the Hamiltonian H in (|J) to ip(x) in ( |l0| ) and 
divide of the exponential we obtain an operator h acting 
on the polynomial Pj_i(x); h has the form 

h = — j^t + (2ix 2 + 2ax + 2ib)-$- 

dx z ax 

-\2i(J - l)x - b 2 - a}. (17) 

Hence, in terms of the generators of the Lie algebra, we 
have 

h = -(J~) 2 + 2iJ+ + 2a J° + 2ibJ- +b 2 + a J. (18) 

This algebraic structure possesses VT symmetry and has 
real eigenvalues. 

We thank the U.S. Department of Energy for financial 
support. 



FIG. 1. The spectrum for the QES Hamiltonian (9) plot- 
ted as a function of b for the case J = 3 and a = 0. For 
b > | (corresponding to the critical value A'criticai = 3) the 
QES eigenvalues are real and are the three lowest eigenvalues 
of the spectrum. When b goes below | , two of the QES eigen- 
values become complex and the third moves into the midst of 
the non-QES spectrum. We believe that the non-QES spec- 
trum is entirely real throughout the (a, b) plane. 



TABLE I. Sequence of critical values for identical and 



-^critical • 


J 


-^critical 


F critical 


J 


-^critical 


^critical 


2 


0.0 


0.0 


16 


25.0526 


-61.3470 


3 


3.0 


-2.0 


17 


26.3475 


-67.3089 


4 


5.47086 


-4.71894 


18 


27.6149 


-73.4116 


5 


7.65570 


-7.93982 


19 


28.8569 


-79.6490 


6 


9.65184 


-11.5572 


20 


30.0754 


-86.0158 


7 


11.5104 


-15.5070 


21 


31.2721 


-92.5072 


8 


13.2625 


-19.7459 


22 


32.4485 


-99.1187 


9 


14.9287 


-24.2419 


23 


33.6058 


-105.846 


10 


16.5235 


-28.9706 


24 


34.7453 


-112.686 


11 


18.0576 


-33.9126 


25 


35.8679 


-119.635 


12 


19.5392 


-39.0521 


26 


36.9747 


-126.689 


13 


20.9747 


-44.3758 


27 


38.0665 


-133.846 


14 


22.3695 


-49.8725 


28 


39.1439 


-141.103 


15 


23.7276 


-55.5323 


29 


40.2078 


-148.458 
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